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L-Ascorbic  acid  (AA),  also  known  as  vitamin  C,  is  an  environmentally-benign  and  biologically-friendly 
compound  that  can  be  used  as  an  alternative  fuel  for  direct  oxidation  fuel  cells.  While  direct  ascorbic  acid 
fuel  cells  (DAAFCs)  have  been  studied  experimentally,  modelling  and  simulation  of  these  devices  have 
been  overlooked.  In  this  work,  we  develop  a  mathematical  model  to  describe  a  DAAFC  and  validate  it 
with  experimental  data.  The  model  is  formulated  by  integrating  the  mass  and  charge  balances,  and  model 
parameters  are  estimated  by  best-fitting  to  experimental  data  of  current- voltage  curves.  By  comparing  the 
transient  voltage  curves  predicted  by  dynamic  simulation  and  experiments,  the  model  is  further  validated. 
Various  parameters  that  affect  the  power  generation  are  studied  by  simulation.  The  cathodic  reaction  is 
found  to  be  the  most  significant  determinant  of  power  generation,  followed  by  fuel  feed  concentration  and 
the  mass-transfer  coefficient  of  ascorbic  acid.  These  studies  also  reveal  that  the  power  density  steadily 
increases  with  respect  to  the  fuel  feed  concentration.  The  results  may  guide  future  development  and 
operation  of  a  more  efficient  DAAFC. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Among  the  various  types  of  fuel  cell,  systems  operating  at  low 
temperature  offer  a  number  of  advantages  that  include  shorter 
start-up  time,  reduced  expense  on  thermal  insulation,  and  security 
[  1,2  ].  The  most  prominent  representative  of  this  group  of  fuel  cells  is 
the  direct  methanol  fuel  cell  (DMFC),  which  has  received  increasing 
attention  recently  due  to  its  easy  storage  of  the  high-energy  den¬ 
sity  liquid  fuel  and  the  simple  reactor  design  that  does  not  require 
fuel  reforming,  as  is  the  case  for  classical  hydrogen  fuel  cells  [3,4]. 
DMFCs  usually  are  adapted  polymer  electrolyte  fuel  cells  (PEMFCs) 
and  operate  with  aqueous  solutions  of  methanol.  Alternative  power 
sources  such  as  DMFCs  and  other  direct-type  PEMFCs  are  expected 
to  be  used  in  portable  electric  devices  for  information  technology 
and  medical  equipment.  However,  the  current  fuels  used  in  direct- 
type  PEMFCs,  namely,  methanol,  ethanol,  propanol  and  ethylene 
glycol,  are  often  toxic  and  this  limits  their  applications  in  medical 
devices,  where  safety  is  a  major  concern. 

Direct  L-ascorbic  acid  (AA)  fuel  cells  (DAAFCs)  may  be  a  solution 
for  this  problem.  L-Ascorbic  acid,  also  known  as  vitamin  C,  is  an 
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environmentally  benign  and  biologically  friendly  compound.  The 
electrochemical  oxidation  of  ascorbic  acid  has  been  the  subject  of 
many  studies  on  various  types  of  electrode  such  as  noble  metals 
and  glassy  carbon  [5-8].  In  addition,  DAAFCs  offer  several  advan¬ 
tages  over  DMFCs  [9].  For  instance,  the  released  product  of  DAAFCs 
is  non-toxic,  and  inexpensive  catalysts  (e.g.,  carbon  black)  can  be 
used  for  oxidation  at  atmospheric  pressure  and  room  temperature. 
While  the  large  crossover  of  methanol  significantly  decreases  the 
cell  voltage  of  the  DMFC  [10],  this  problem  is  not  serious  in  DAAFC 
as  experiments  show  that  there  is  very  little  crossover  of  L-ascorbic 
acid  through  a  Nafion  membrane  of  a  certain  thickness  [9]. 

While  the  characteristics  of  DAAFCs  have  been  studied  by  exper¬ 
iments  and  the  maximum  power  output  was  found  to  be  about 
one-third  that  of  DMFCs  [11],  the  performance  of  the  DAAFCs  can 
be  further  investigated  by  simulation.  The  performance  of  a  DAAFC, 
like  its  counterpart  DMFC,  is  influenced  by  many  parameters. 
Investigating  each  parameter  independently  and/or  all  possible 
combinations  of  parameters  using  laboratory  experiments  would 
be  costly  and  time-consuming.  Numerical  modelling  of  a  fuel  cell 
system  is  a  valuable  research  tool  for  investigating  system  param¬ 
eters  that  would  reduce  the  amount  of  time  and  money  required, 
because  of  the  ease  with  which  models  can  be  modified  to  simulate 
various  configurations  and  operating  conditions.  While  modelling 
and  simulation  has  demonstrated  its  importance  in  developing 
DMFCs  and  other  chemical  fuel  cell  systems  [3,10,12-19],  it  remains 
overlooked  in  DAAFC  development.  In  order  to  utilize  DAAFCs  in 
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Fig.  1.  Schematic  diagram  of  the  direct  ascorbic  acid  fuel  cell. 


an  effective  way  and  improve  their  overall  performance,  a  math¬ 
ematical  model  is  necessary  so  that  the  system  can  be  analyzed 
and  optimized  by  simulation  for  different  operating  conditions.  The 
capability  of  predicting  transient  dynamics  subject  to  the  changes 
in  the  power  demands  and  exterior  conditions  are  of  high  inter¬ 
est  for  mobile  applications.  Therefore,  the  purpose  of  the  present 
work  is  to  develop  a  mathematical  model  for  DAAFCs  which  allows 
not  only  studying  steady-state  behaviour  but  also  investigating 
dynamic  responses. 

The  proposed  model  is  implemented  in  MATLAB  7.0  [20].  The 
model  is  developed  on  the  basis  of  mass  and  charge  balances,  and 
it  describes  the  electrochemical  kinetics  in  the  catalyst  layers  and 
the  mass  transfer  between  the  bulk  solutions  and  catalyst  layers. 
Steady-state  data  from  experimentally  measured  current-voltage 
curves  are  used  for  fitting  the  model  parameters  and  validating 
the  model.  The  dynamic  behaviour  of  step  changes  and  periodic 
changes  in  fuel  flow  rate  and  concentration  are  compared  for  exper¬ 
imental  and  simulated  results.  The  sensitivity  analysis  of  various 
model  parameters  is  carried  out  by  simulation  and  the  most  sig¬ 
nificant  factors  that  affect  power  generation  are  discussed.  These 
results  promise  to  be  useful  for  developing  or  operating  more  effi¬ 
cient  DAAFCs. 

2.  Methods  of  approach 

2.1.  Experiments 

A  schematic  diagram  of  the  experiment  set-up  of  the  DAAFC 
is  presented  in  Fig.  1.  The  fuel  cell  was  fitted  with  a  mem¬ 
brane  electrode  assembly  (MEA)  that  consisted  of  a  NAFION 


Table  1 

Experimental  conditions  and  constants 


Parameter 

Description 

Units 

Value 

U° 

Open  circuit  voltage 

V 

0.592 

T 

Cell  operating  temperature 

°C 

25 

po2 

02  pressure 

atm 

1 

As 

Cross-sectional  electrode  area 

m2 

0 

X 

0 

km 

Conductivity  of  membrane 

£2_1  m-1 

17 

dm 

Thickness  of  the  membrane 

m 

2  x  10-4 

Va 

Anode  compartment  volume 

m3 

2.0  x  10-6 

F 

Faraday  constant 

Cmol-1 

96485.4 

R 

Gas  constant 

Jmol-1  K-1 

8.3144 

yCL 

Void  volume  of  anode  catalyst  layer 

m3 

8.0  x  10-9 

Ca 

Anode  capacitance 

Fm~2 

400 

Cc 

Cathode  capacitance 

Firr2 

500 

117™  PEM  in  the  middle,  anodic  and  the  cathodic  catalyst  lay¬ 
ers  on  either  side,  diffusion  layers,  and  the  reactant  distribution 
channels.  Fuel  cell  experiments  were  performed  using  a  sin¬ 
gle  cell  at  atmospheric  pressure  and  room  temperature  (25  °C). 
The  parameters  of  the  experimental  set-up  are  summarized  in 
Table  1. 

L-Ascorbic  acid  was  supplied  as  an  aqueous  solution  to  the 
anode,  and  oxygen  saturated  with  water  was  fed  at  100  ml  min-1  to 
the  cathode  as  the  oxidant.  In  order  to  estimate  the  model  param¬ 
eters  with  better  accuracy,  the  experiments  were  carried  out  by 
combining  different  levels  of  AA  concentrations  and  flow  rates. 
Specifically,  the  concentration  was  selected  in  the  range  of  0.1-1  M 
with  the  upper  limit  close  to  the  AA  solubility  and  the  flow  rate 
in  the  range  of  1-8  ml  min-1,  which  covers  the  operating  ranges 
of  the  DAAFC.  The  experiments  were  then  performed  with  combi¬ 
nations  of  high,  middle  and  low  levels  of  concentration  and  high, 
middle  and  low  flow  rates.  As  a  result,  nine  runs  of  experiments, 
as  listed  in  Table  2,  were  conducted  and  a  total  of  nine  polariza¬ 
tion  (current- voltage)  curves  were  measured.  All  experiments  were 
carried  out  under  galvanostatic  conditions,  i.e.,  the  electrical  cell 
current  was  fixed  as  an  input  variable  and  the  resulting  cell  voltage 
was  measured  as  an  output  variable. 

To  validate  further  the  developed  DAAFC  model  and  gain  better 
insight  into  the  fuel  cell  performance,  dynamic  experiments  were 
performed.  Similar  to  the  study  of  DMFC  dynamics  [3,10],  the  inves¬ 
tigation  of  the  dynamic  behaviour  of  the  DAAFC  was  undertaken  by 
holding  certain  variables  constant  while  observing  other  variables, 
e.g.,  the  cell  voltage  or  concentrations  of  reactants.  Specifically,  the 
excitation  functions,  i.e.,  step  function,  and  periodic  function  of  con¬ 
trolled  variables,  namely,  AA  feed  concentration  and  flow  rate,  were 
applied,  and  the  resultant  responses  (e.g.,  the  resultant  variation  of 
the  cell  voltage  with  time)  were  measured  and  used  to  compare 
with  the  model  predictions. 


Table  2 

Variations  of  experimental  conditions  and  the  resultant  mass  transfer  coefficients 


Data  set  # 

Concentration  of 
ascorbic  acid  (AA) 

(mol  L-1 ) 

AA  feed  flow  rate  (ml  min-1 ) 

Limiting  current 
density  readings 
(mA  cm-2) 

Mass  transfer 
coefficient  I^s 
(10-5  ms-1) 

Average  voltage 
difference  between 
the  experiment  and 
simulation  (V) 

1 

0.5 

4 

95.1 

1.1566 

0.008 

2 

0.5 

1 

63.3 

1.0820 

0.011 

3 

0.5 

8 

108 

1.2219 

0.005 

4 

0.1 

4 

30.3 

2.0539 

0.012 

5 

0.1 

8 

40.0 

2.4544 

0.008 

6 

1 

1 

85.6 

0.6045 

0.010 

7 

1 

4 

122 

0.6985 

0.007 

8 

1 

8 

133 

0.7268 

0.014 

9 

0.2 

8 

90 

2.8263 

0.008 
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2.2.  Assumptions  of  modelling 

As  described  in  Section  2.1,  the  basic  concepts  of  modelling  a 
single-cell  DAAFC  encompasses  a  PEM  in  the  middle,  the  anodic 
and  the  cathodic  reaction  zones  (catalyst  layers)  on  both  sides,  and 
anode  and  cathode  compartments  (diffusion  layers  and  the  reactant 
distribution  channels).  Since  the  design  configuration  and  working 
principle  of  the  DAAFC  are  similar  to  that  of  a  DMFC  [10],  similar 
assumptions  were  made  as  follows: 


In  the  cathode,  the  reduction  of  dissolved  oxygen  is  described  as 
(R2).  Our  preliminary  study  found  that  the  reversed  reaction  was 
negligible.  Thus,  the  reaction  rate  of  (R2)  is  formulated  as  Eq.  (2), 
where  P02  is  the  oxygen  partial  pressure  in  the  oxygen  saturated 
with  water  at  25  °C  and  fed  to  the  cathode  compartment,  rjc  is  the 
cathode  overpotential,  designates  the  rate  constant  of  the  cath¬ 
ode  reaction  at  standard  conditions,  and  is  the  charge-transfer 
coefficient  of  cathode.  It  is  noted  that  r2  is  defined  in  terms  of 
electron  generation,  just  like  r\  in  Eq.  (1). 


•  The  anode  compartment  is  treated  as  a  continuous  stirred  tank 
reactor  (CSTR). 

•  Ohmic  drops  in  current-collectors  and  electric  connections  are 
negligible. 

•  The  fuel  cell  operates  isothermally. 

•  Oxygen  is  fed  in  excess,  i.e.,  oxygen  conversion  in  the  cathode 
compartment  is  negligible,  and  therefore  no  oxygen  mass  balance 
is  required. 

•  Oxygen  does  not  diffuse  into  the  PEM. 

•  Water  concentration  is  constant  (excess  component  in  liquid  mix¬ 
ture). 

•  Mass-transport  resistances  in  the  catalyst  layers  are  negligible. 

•  Mass-transport  coefficients  of  ascorbic  acid  (AA)  and  its  product 
dehydroascorbic  acid  (DHA)  in  the  anode  diffusion  layer  are  equal. 

•  The  mixtures  in  the  anode  compartment  and  in  the  anode  diffu¬ 
sion  layer  are  treated  as  pure  liquids. 

2.3.  Kinetics ,  mass  and  charge  balances 

Experiments  suggested  that  ascorbic  acid  is  oxidized  in  the 

anode  by  a  two-electron  transfer  reaction,  i.e.,  the  same  as  its 

metabolic  conversion  [  1 1,21,22 ],  which  is  described  as  (R1 ).  The  rate 

expression  of  the  reaction  takes  the  Butler-Volmer  form,  as  shown 


by: 

H 

H 

HCOH 

HCOH 

1 

HCOH 

HCOH 

^V°V=o 

>  w 

-r\  _i_  o  u+  _i_  o  q- 

\_r 

-U  +  2  rl  +  2  e 

HO  OH 

n 

o  o 

(AA) 

(DHA) 

(Rl) 

r,  =k°exp(|^a) 

c^-^exp(^. 

\cCL 
j  ^DHA 

(1) 

0.502  +  2H+  +  2e_  H20 


(R2) 


r2  =  -k°(P02)1/2exp 


06—1  )F 

RT 


he 


(2) 


In  a  similar  manner  to  the  work  of  Sundmacher  et  al.  [3],  where 
the  dynamics  of  a  direct  methanol  fuel  cell  were  studied,  the  gov¬ 
erning  equations  associated  with  the  mass  balances  of  AA  and  DHA 
can  be  formulated  in  Eqs.  (3)-(8)  with  the  assumption  that  the 
anode  is  a  continuous  stirred  tank  reactor  (CSTR)  and  functions  as 
a  storage  module.  In  these  equations,  the  subscripts  ‘a’  indicates 
the  anode  compartment,  the  superscript  ‘in’  the  feed  flow,  and  Va, 
V^L,  Q.  and  As  are  the  volumes  of  the  anodic  compartment  and  cat¬ 
alyst  layer,  the  AA  feed  flow  rate,  and  the  cross-section  area  of  the 
membrane,  respectively. 


dCAA  _  Q  rrin 

dt  “ValAA 


Caa) 


/clsas 

”vr" 


(Caa  -  Caa) 


(3) 


Cl^HA  =  ^(CDHA  “  CDHA)  -  — yT~(CDHA  “  C^A)  (4) 

The  second  terms  on  the  right-hand  side  of  Eqs.  (3)  and  (4), 
respectively,  represents  the  mass  fluxes  of  AA  and  DHA  in  the  anode 
diffusion  layer  where  no  storage  is  assumed.  It  is  noted  that  the 
mass-transfer  coefficient  I^s  is  taken  to  be  equal  for  both  AA  and 
DHA.  In  the  catalyst  layer,  the  mass  balances  for  AA  and  DHA  in  the 
catalyst  layer  are  given  by  Eqs.  (5)  and  (6). 


deft  J<lsAS  ~  As 

dt  ~~  yCL  ^AA  Caa)  yCl/1 


(5) 


dcCL 

U^DHA 

dt 


JClsAs  _  CL  ,  As 

yCL  (^DHA  -  ft) ha)  +  yCL  rl 


(6) 


The  charge  balances  at  the  anode  and  cathode  are  given  as  below 
where  icen  is  the  cell  current  density,  Ca  and  Cc  are  charge  double¬ 
layer  capacitances  of  the  anode  and  the  cathode,  respectively,  which 
are  measured  by  experiment. 


where  Caa  and  CDHa  are  the  concentrations  of  AA  and  DHA,  respec¬ 
tively;  the  superscript  ‘CL’  indicates  the  catalyst  layer;  rja  is  the 
anode  overpotential.  In  addition,  k°  indicates  the  rate  constant  of 
the  anode  reaction  at  standard  conditions,  K\  is  the  equilibrium 
constant,  a  is  the  charge-transfer  coefficient  of  the  anode,  F  is  the 
Faraday  constant,  R  is  the  gas  constant,  and  T  is  the  operating  tem¬ 
perature  of  the  cell.  Notably,  a  reversible  reaction  is  assumed  for 
the  AA  oxidation. 


^  =  T(icel,-2Fr1)  (7) 

^£  =  T(-icell-2Fr2)  (8) 

As  mentioned  previously,  the  ohmic  drops  in  the  current- 
collectors  and  electric  connections  are  negligible.  Thus  the  cell 
voltage  Uceii  can  be  calculated  as  below,  where  U°  is  the  open-circuit 


Table  3 

Resultant  parameter  values 


Symbol 

Description 

Unit 

Value 

*? 

Forward  rate  constant  of  AA  oxidation  (Rl ) 

molm-2  s-1 

2.91  x  1CT6 

Kl 

The  ratio  of  forward  rate  to  backward  rate  of  AA  oxidation  (Rl ) 

- 

3.11  x  10-2 

k2 

Standard  rate  constant  of  02  reduction  (R2) 

molm-2  s-1 

4.47  x  10~6 

a 

Charge  transfer  coefficient  of  anodic  reaction 

- 

0.410 

P 

Charge  transfer  coefficient  of  cathodic  reaction 

- 

0.467 

98 
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potential,  and  km  and  dm  are  the  conductivity  and  thickness  of  the 
membrane,  respectively. 

Uceii  =  [/°-??a  +  >7c-  (^)l'cell  O) 

2.4.  Parameter  estimation 

In  order  to  make  use  of  the  developed  model  (Eqs.  (l)-(9))  to 
predict  the  DAAFC  performance,  the  associated  parameters  or  con¬ 
stants  have  to  be  determined  a  priori.  Some  of  the  parameters  have 
been  measured  by  experiments,  as  shown  in  Table  1,  while  oth¬ 
ers  are  determined  by  the  mathematical  methods  described  in  this 
section. 


2.4A.  Estimation  of  mass-transfer  coefficients 

The  effective  mass-transfer  coefficient  I<ls  is  determined  as 
follows  with  the  measured  limiting  current  density.  In  the  steady- 
state,  Eqs.  (3),  (5)  and  (7)  are  all  zero,  and  the  AA  concentration  in 
anode  catalyst  layer  is  also  zero  at  limiting  current.  Combining  Eqs. 


(3),  (5)  and  (7)  yields 

KLS  =  - Ohm -  (10) 

2FQC^-A%m 

As  can  be  seen  from  the  above  equation,  in  addition  to  the  lim¬ 
iting  current  density,  the  mass-transfer  coefficient  is  also  affected 
by  both  the  AA  concentration  and  the  flow  rate  of  the  feed.  The 
measured  limiting  current  densities  are  listed  in  Table  2  (column 
4),  and  the  resultant  mass-transfer  coefficients  of  the  nine  runs  of 
steady-state  experiments  are  evaluated  and  presented  in  Table  2 
(column  5). 

2.4.2.  Estimation  of  model  parameters 

As  stated  in  Ref.  [23],  where  a  DMFC  model  is  studied,  the  model 
parameters  are  difficult  to  determine  directly  from  experiments. 
These  parameters,  as  listed  in  Table  3,  are  /<i,  k®,  a  and  /3.  The 
model  parameters  are  estimated  in  the  following  manner  [24]. 
Specifically,  the  model  parameters  are  estimated  by  minimizing  the 
squared  differences  of  the  voltages  over  six  sets  of  data,  namely  the 
quantity  Fobj  defined  in  Eq.  (11 )  where  N  is  the  number  of  data  sets 
(N= 6),  and  M;  designates  the  number  of  data  points  of  the  ith  data 


Run-2:  C=0.5M,  Q=1  ml  min'1 


Run-3:  C=0.5M,  Q=8  ml  min'1 


0  100  200  300  400  500  600  700  800 

Current  Density  (A  m'2) 


Run-8:  C=1M,  Q=8  ml  min'1  Run-9:  C  =  0.2M,  Q  =  8ml  min'1 


Fig.  2.  Comparisons  of  the  experimental  measurements  and  simulated  results  from  best  fitting  the  model  using  data  sets  #2,  #3,  #5,  #6,  #8,  and  #9. 
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Table  4 

Average  percentage  changes  in  power  output  due  to  respective  change  of  each 
parameter  by  the  factor  of  0.9  and  1.1 


Ratio  p 

Jds 

k°2 

a 

p02 

*? 

I< 1 

Q 

dm 

0.9 

11.6 

-5.6 

-3.6 

-2.5 

-2.6 

-1.3 

-1.2 

-0.9 

-0.6 

0.3 

1.1 

-13.8 

4.0 

2.6 

2.3 

1.9 

1.1 

1.0 

0.8 

0.5 

-0.3 

set. 


f*-E  S  m, 


1=1  \j=  1 


(11) 


The  simulated  cell  voltage  U9J1  is  obtained  from  Eq.  ( 9 )  by  solving 
the  system  of  model  Eqs.  (l)-(8)  at  the  steady  state,  i.e.,  the  alge¬ 
braic  equations  obtained  by  setting  all  the  derivatives  to  zeros.  It  is 
noted  that  six  sets  of  experimental  measurements,  i.e.,  sets  #2,  #3, 
#5,  #6,  #8,  #9,  are  used  to  fit  the  model  and  estimate  the  parame¬ 
ter  values,  while  the  remaining  three  sets,  namely  #1,  #4,  #7,  serve 
for  validation  of  the  developed  model.  The  resultant  parameters 
are  listed  in  Table  3.  Comparisons  of  the  experimental  measure¬ 
ments  and  simulated  voltages  for  the  six  fitting  sets  are  shown  in 
Fig.  2.  As  can  be  seen  from  the  plots,  good  agreement  is  observed 
between  the  experimental  and  simulated  results  generated  from 
the  best-fitting  model.  Furthermore,  the  model  validation  is  also 
successful,  as  shown  in  the  three  plots  of  Fig.  3.  The  average  dif¬ 
ferences  between  the  experiments  and  model  predictions  are  in 
the  range  of  0.008-0.014  (V),  as  shown  in  column  6  of  Table  2.  It 
is  noted  that  the  differences  are  generally  below  the  acceptable 
measurement  error  of  10%  that  is  equivalent  to  the  minimum  mea¬ 
surement  error  of  0.01  V.  Since  the  three  sets  of  validation  data 
are  not  involved  in  fitting  the  model,  it  may  be  concluded  that  the 
resulting  model  demonstrates  quite  good  predictive  capability. 


3.  Results  and  discussion 

3.1  Parameter  sensitivity  analysis 

A  primary  sensitivity  analysis  of  the  operating,  design  and  model 
parameters  was  investigated  by  varying  each  parameter  in  turn, 
while  leaving  the  other  parameters  unchanged.  The  average  per¬ 
centage  changes  in  the  power  output  over  a  range  of  cell  current 
density  of  50-700  (A  nrr  2 )  are  listed  in  Table  4,  where  each  param¬ 
eter  is  varied  by  10%.  In  other  words,  the  parameter  is  respectively 
multiplied  by  a  ratio  of  0.9  and  1.1  to  its  original  value  at  the  base 
case,  i.e.,  the  experimental  condition.  As  shown  in  Table  4,  when  the 
parameter  is  reduced  by  a  factor  of  0.9,  the  electron-transfer  coef¬ 
ficient  of  cathode,  /3,  results  in  the  largest  boost  in  power  density 
by  an  average  of  11.6%  over  the  entire  range  of  cell  current  den¬ 
sity.  The  next  largest  change  is  the  AA  feed  concentration  which 
results  in  a  decrease  of  5.6%.  The  mass  transfer  coefficient  I^s  ranks 
the  third  significant  factor  which  causes  a  3.6%  decrease  in  power 
output.  The  order  of  significance  of  the  parameters,  as  shown  by 
the  absolute  values  of  the  percentages,  is  listed  in  Table  4  with  the 
highest  on  the  left  and  the  lowest  on  the  right.  It  is  noted  that  in  the 
case  of  a  reduction  of  parameters,  only  and  membrane  thickness 
lead  to  an  increase  in  power  output,  while  others  decrease  it.  On 
the  other  hand,  the  order  of  parameter  significance  remains  almost 
unchanged  as  the  parameters  are  increased  by  a  factor  of  1.1,  but 
opposite  changes  in  power  generation  are  observed.  In  both  cases, 
the  membrane  thickness  is  the  least  significant.  A  visual  exam¬ 
ination  of  the  resultant  power  output  due  to  changes  in  all  the 
parameters  can  be  found  in  Fig.  4.  As  can  be  seen,  higher  current 
densities  results  in  larger  deviations  from  the  base  case. 


Validation  Run-1:  C=0.5M,  Q=4  ml  min'1 


Validation  Run-4:  C=0.1M,  Q=4  ml  min'1 


Fig.  3.  Model  validation  by  comparing  the  experimental  measurements  and  the 
model  predictions  using  data  sets  #1,  #4,  #7. 


The  sensitivity  of  the  five  model  parameters  listed  in  Table  3  was 
carried  out  by  the  local  relative  sensitivity  analysis  method  [25],  to 
evaluate  the  ratio  of  changes  in  the  computed  power  density  to 
changes  in  the  parameter  values.  The  following  equation  was  used 
for  the  five  model  parameters. 


py^+s^-p^xj)  x, 

8xj  P{t,Xj)  J 


(12) 


where  Tj  is  time-dependent  sensitivity  of  the  jth  parameters;  Xj 
is  the  value  of  the  jth  parameter;  8xj  is  the  change  in  xf,  P  is  the 
power  density.  In  the  present  study,  8xj  =  5 %Xj.  A  dynamic  simu¬ 
lation  of  step  changes  in  AA  feed  flow  rate  from  8  to  4  and  then 
back  to  8  ml  min-1  was  conducted  to  examine  the  sensitivities  of 
the  five  parameters.  The  higher  the  sensitivity  of  the  parameter, 
the  farther  the  sensitivity  curve  locates,  apart  from  zero.  As  shown 
in  the  results  (Fig.  5),  f3  sensitivity  is  the  farthest  and  stands  out 
from  all  of  the  rest.  Again,  this  indicates  that  is  the  most  sensitive 
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Fig.  4.  The  power  density  resulted  from  (a)  decrease  in  parameters  by  the  factor  of 
0.9  and  (b)  increase  in  parameters  by  the  factor  of  1.1. 


parameter.  The  order  of  the  sensitivity  of  the  remaining  four  param¬ 
eters  (from  highest  to  the  lowest)  is  a  >  >  /<?  >  K\.  Overall,  the 

cathodic  reaction  is  the  most  significant  factor  limiting  the  cell  per¬ 
formance,  followed  by  the  AA  feed  concentration  and  mass-transfer 
coefficient. 

3.2.  Steady-state  analysis 

With  the  resulting  parameters  estimated  in  Section  2.4,  the  com¬ 
ponent  concentrations  and  overpotentials  of  the  anode  and  cathode 
were  evaluated  with  respect  to  the  changes  in  cell  current  density. 
Since  the  resultant  trends  of  each  run  of  simulation  are  similar,  only 
the  results  from  run  #1  are  presented  in  this  paper.  As  shown  in 
Fig.  6a,  the  AA  concentrations  in  both  the  anode  compartment  and 
the  catalyst  layer  linearly  decrease  with  an  increase  in  cell  current 
density,  whereas  the  DHA  concentrations  increase  linearly.  This 
can  be  explained  by  the  fact  that  higher  current  density  leads  to 
a  larger  reaction  rate,  i.e.,  faster  consumption  of  the  reactant  AA 
and  production  of  DHA.  Furthermore,  the  AA  concentration  in  the 
compartment  is  higher  than  that  in  the  catalyst  layer  due  to  the 
mass-transfer  limitation.  For  the  same  reason,  the  DHA  generated 
in  the  catalyst  layer  demonstrates  opposite  trends  to  those  of  AA 
concentration.  In  Fig.  6b,  it  is  shown  that  the  anode  overpotential 
increases  while  the  cathode  overpotential  decreases  with  respect 
to  an  increase  of  current  density.  These  results  explain  the  decrease 
of  cell  voltages  (Fig.  2). 

As  mentioned  in  the  preceding  section,  the  fuel  (AA)  feed 
concentration  is  the  second  significant  factor  affecting  power  gen¬ 
eration.  We  therefore  studied  further  the  concentration  effect  by 
simulation.  In  particular,  the  simulation  was  performed  with  the 
AA  feed  concentration  varying  from  0.1  to  1.9  M  while  the  flow 
rate  remained  unchanged  at  4  ml  min-1  and  other  parameters  also 
remained  unchanged.  It  is  observed  in  Fig.  7a  that  the  power  density 
increases  with  increase  in  fuel  feed  concentration,  and  the  attain¬ 
able  current  density  range  also  increases.  The  optimum  current 
density  that  results  in  maximum  power  density  with  respect  to 
the  fuel  feed  concentration  is  plotted  in  Fig.  7b. 

3.3.  Dynamic  analysis 
3.3 A.  Step  changes 

In  this  work,  we  have  studied  the  dynamic  behaviour  of  the 
DAAFC  by  simulation  and  then  compared  the  results  with  those 
from  experiments.  Specially,  step  changes  in  AA  feed  concentra¬ 
tion  and  flow  rate  are  carried  out  respectively.  For  the  former  case, 
a  step  decrease  in  AA  feed  concentration  occurs  from  0.5  to  0.2  M  at 
time  =  100  s  with  a  flow  rate  of  AA  feed  at  8  ml  min-1  and  the  cur¬ 
rent  density  maintained  at  100,  200  and  300 Am-2,  respectively. 
Operation  prior  to  the  step  change  (i.e.,  C=0.5M)  is  equivalent  to 
the  steady-state  run  #3.  The  simulation  and  experiment  results  are 
given  in  Fig.  8a.  In  general,  the  model  is  able  to  simulate  the  experi¬ 
ment  although  deviations  are  found.  For  instance,  a  relatively  larger 
discrepancy  is  observed  at  current  density  of  100  Am-2.  It  is  noted 
that  the  discrepancy  results  from  fitting  the  model  parameters  with 
all  the  six  sets  of  experimental  readings.  As  can  be  seen  from  the  fit¬ 
ting  result  of  run  #3  (C  =  0.5  M  and  flow  rate  =  8  ml  min-1 )  in  Fig.  8b 
that  is  reproduced  from  Fig.  2,  point  A  which  corresponds  to  A  in 
Fig.  8a  is  not  well  fitted.  In  fact,  point  A  appears  somewhat  anoma¬ 
lous  in  the  set  of  experimental  readings  that,  otherwise,  are  quite 
smooth.  Therefore,  the  less  consistent  fitting  of  experiment  and 
simulation  may  lead  to  the  observed  discrepancy  in  the  step-change 
curve.  By  contrast,  point  B  (200  Am-2)  and  point  C  (300  Am-2)  are 
well  fitted  in  Fig.  8b,  thus  good  agreements  are  observed  in  the  cor¬ 
responding  step-change  curves  (Fig.  8a).  It  is  noted  that  there  is  a 
valley  in  each  voltage  curve,  which  is  due  to  the  practical  switch  of 
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Fig.  6.  Steady-state  simulation  result  for  run  #1.  (a)  AA  and  DHA  concentrations  in  anode  compartment  and  catalyst  layer  and  (b)  anode  and  cathode  overpotentials. 


Fig.  7.  (a)  Simulated  power  density  curves  at  various  AA  feed  concentrations  and  (b)  maximum  power  density  with  respect  to  different  AA  feed  concentrations. 


the  feed  flow  concentration,  i.e.,  in  the  experiment,  the  pump  was 
stopped  for  5-6  s  for  each  switch. 

In  the  latter  case,  a  step  change  in  AA  feed  flow  rate  from  8 
to  1  ml  min-1  was  made  at  time  =  50  s,  with  a  constant  AA  feed 
concentration  of  0.5 M  and  a  current  density  200  and  400 Am-2, 
respectively.  This  is  equivalent  to  switching  from  steady-state  run 
#3  to  run  #2  (see  Fig.  2).  The  simulation  and  experimental  results 
are  shown  in  Fig.  9a.  A  minor  discrepancy,  as  shown  by  point  D, 
can  be  found  for  a  current  density  of 400  A  m-2  before  the  flow  rate 


is  switched.  This  point  is  equivalent  to  point  D  in  Fig.  8b  which  is 
the  same  as  the  steady-state  result  shown  in  Fig.  2  (run  #3)  where 
fairly  good  agreement  is  observed.  Further  examination  reveals  that 
the  discrepancy  is  caused  by  the  inconsistent  results  of  the  two 
methods  of  voltage  measurement.  In  other  words,  all  the  steady- 
state  I-V  data  (i.e.,  run  #1  to  #9)  are  measured  by  a  current  scan 
and  all  the  dynamic  I-V  data  by  a  voltage  scan.  In  this  work,  the 
model  parameters  are  estimated  by  best-fitting  all  the  steady-state 
measurements,  and  dynamic  data  serve  as  further  validation  of  the 


Current  Density  (A  m'2) 


Fig.  8.  (a)  Comparison  of  cell  voltage  of  the  dynamic  simulation  (fine  lines)  and  experimental  (bold  lines)  results  for  the  step  changes  of  AA  feed  concentration  from  0.5  to 
0.2  M  with  constant  AA  feed  flow  rate  of  8  ml  min-1 .  (b)  Steady-state  simulation  and  experiment  results  at  AA  concentration  =  0.5  M  and  flow  rate  =  8  ml  min-1  (run  #3). 
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Fig.  9.  (a)  Comparison  of  cell  voltage  of  the  dynamic  simulation  (fine  lines)  and  experimental  (bold  lines)  for  the  step  changes  of  AA  feed  flow  rate  from  8  to  1  ml  min-1  with 
constant  AA  feed  concentration  of  0.5  M.  (b)  Steady-state  simulation  and  experiment  results  at  AA  concentration  =  0.5  M  and  flow  rate  =  1  ml  min-1  (run  #2). 


developed  model.  This  may  also  explain  the  discrepancy  occurring 
at  400  A  m-2  and  1  ml  min-1  (corresponding  to  run  #2).  As  seen  at 
point  B  in  Fig.  9b,  good  agreement  is  obtained  between  the  simu¬ 
lation  and  the  steady-state  experiment  as  measured  by  a  current 
scan,  but  measurements  by  a  voltage  scan  for  the  dynamic  exper¬ 
iment  oscillate  and  the  average  is  slightly  lower  than  that  of  the 
steady  state.  During  200  Am-2  and  1  ml  min-1  (also  correspond¬ 
ing  to  run  #2),  the  deviation  may  result  from  the  abovementioned 
issue  or  it  may  inherit  from  the  mismatch  during  fitting  the  model 
parameters  at  steady  state  (point  A  in  Fig.  9b). 

3.3.2.  Periodic  changes 

In  the  experiment  of  periodic  dynamics,  the  delivery  pump  was 
switched  off  at  the  18  s  point  and  it  took  5-6  s  to  stop  the  action 
completely.  The  solution  was  changed  from  0.5  M  AA  to  FI20  and 
the  pump  was  switched  on  immediately.  The  water  feed  started  at 
around  24  s.  It  is  estimated  that  it  took  about  2  s  to  change  the  fuel 
solution  from  0.5  M  AA  to  water.  The  fuel  solution  was  changed  at 
24  s  intervals  by  repeating  the  above  operation.  Since  the  feed  solu¬ 
tion  in  the  tank  approaches  the  inlet  of  the  fuel  cell  via  three  tubes 
and  a  pump  with  a  relatively  large  volume,  the  switch  in  fuel  solu¬ 
tion  concentration  (from  0.5  M  to  H20)  cannot  immediately  affect 
the  fuel  cell  operation.  In  fact,  the  tubes  and  pump  acted  like  a 
buffer.  The  original  AA  concentration  of  the  feed  can  be  regarded  as 
unchanged  (0.5  M)  after  the  first  few  switches  until  the  solution  in 
the  tubes  is  replenished.  Thus  the  0.5-M  solution  sustains  the  fuel 
cell  operation  for  the  first  few  cycles  of  pump  stopping.  Thereafter, 
the  solution  in  the  tubes  is  diluted.  The  experiments  were  carried 
out  at  current  density  of  100,  200  and  300 Am-2  and  the  AA  flow 
rate  was  maintained  at  8  ml  min-1 . 

By  estimating  the  volumes  of  the  buffer  and  the  cell,  it  was  found 
that  during  the  first  four  cycles  AA  was  fed  at  0.5  M.  Because  of 
the  presence  of  a  buffer  to  the  AA  feed  flow,  the  concentration 
was  switched  from  0.5  M  to  an  average  diluted  value,  which  was 
calculated  as  0.44  M.  A  comparison  of  the  simulation  and  experi¬ 
mental  results  is  shown  in  Fig.  10.  The  simulation  tracks  the  periodic 
changes  in  the  voltage  curves,  but  in  a  more  stable  fashion.  A  rela¬ 
tively  large  deviation  is  found  at  a  current  density  of  100  A  m-2  and 
is  thought  to  be  due  to  a  previous  mismatch  in  fitting  the  steady 
state  (point  A  in  Fig.  8b).  It  is  observed  that  some  experiment  mea¬ 
surements  demonstrate  very  sharp  valleys  while  the  simulation 
does  not  generate  such  results.  The  very  low  experimental  volt¬ 
ages  are  probably  caused  by  the  very  low  feed  concentration  when 
water  is  being  pumped.  It  can  be  envisaged  that  the  concentration 


Time  (Sec) 

Fig.  10.  Comparison  of  cell  voltage  of  the  dynamic  simulation  (fine  lines)  and  exper¬ 
imental  (bold  lines)  for  the  periodic  changes  of  AA  feed  flowrate  with  constant  AA 
feed  flow  rate  of  8  ml  min-1 . 

change  prior  to  the  inlet  of  the  cell  is  quite  complicated  and  depends 
not  only  on  the  amount  of  the  water  but  also  on  the  geometry  of 
the  buffer  space.  Since  the  focus  has  been  on  modelling  of  the  fuel 
cell  unit  rather  than  of  the  feeding  facilities,  we  cannot  capture  the 
patterns  in  feed  concentration.  In  this  work,  the  average  feed  con¬ 
centration  has  been  taken  for  the  calculation  and  thus  very  low 
voltages  are  not  predicted. 

4.  Conclusions 

The  present  work  demonstrates  a  modelling  method  to  describe 
direct  ascorbic  acid  fuel  cells  (DAAFCs).  In  general,  the  simula¬ 
tion  results  agree  with  the  experimental  data  at  steady  or  dynamic 
states.  Since  the  present  study  focuses  only  on  the  fuel  cell  itself, 
however,  changes  in  the  feed  flow  concentration  prior  to  approach¬ 
ing  the  inlet  cannot  be  accurately  simulated.  Furthermore,  the 
assumption  of  ideal  mixing  in  the  anode  compartment  may  also 
cause  less  accurate  prediction  because  the  anode  is  a  narrow  and 
relatively  long  channel.  As  a  result,  the  transient  behaviour  of  the 
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simulated  voltages  does  not  completely  match  that  of  experiments 
at  the  points  where  a  step  change  occurs. 

The  parameter  sensitivity  analysis  by  simulation  reveals  that 
the  cathode  electron  transfer  coefficient  is  the  most  significant 
factor  limiting  the  cell  performance,  followed  by  the  AA  feed  con¬ 
centration  C™A  and  mass-transfer  coefficient  I^s.  To  improve  the 
DAAFC  performance  by  varying  /3  and  /CLS,  the  important  issues  are 
to  develop  more  efficient  cathode  materials  and  catalyst  materials. 
This  is,  however,  beyond  the  scope  of  the  present  study  and  should 
be  the  subject  of  further  investigation.  Simulation  shows  that  an 
increase  in  CJ^  will  lead  to  a  boost  of  the  power  generation.  This 
result  may  imply  that  in  order  to  obtain  higher  power  generation, 
the  AA  feed  concentration  should  be  increased  as  high  as  possible. 
The  present  investigation  offers  further  insight  into  how  various 
parameters  affect  the  power  output,  and  may  provide  guidance  for 
more  efficient  DAAFC  design  and  operational  strategies. 
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